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Abstract 

Color superconductivity of QCD at finite density, temperature and flavor asym- 
metry is studied within an approximation which the interaction is modeled upon 
four-fermion interactions. We calculate the thermodynamic potential in the so- 
called NJL-type model and study the phase structure at finite density, temperature 
and flavor asymmetry. We find that a mixed phase appears at sufficiently large 
flavor asymmetry and low temperature. A tricritical point in the T-^B-fJ-i plane 
is also found. 

1 INTRODUCTION 

Quantum Chromodynamics (QCD) describes the interaction of quarks and gluons. Since 
QCD is asymptotically free, the coupling constant becomes small when either tempera- 
ture or the chemical potential, namely the the Fermi momentum, is large. By Cooper 
theorem, Fermi surface has instability when an arbitrarily small attractive interaction 
between quarks presents. For quark-quark scattering, where the quarks are in the an- 
titriplet 3, one gluon exchange interaction is attractive and in the 6, it is repulsive. 
Thus we expect that, in the case of two quark flavors, the quark-quark condensate in 3 
channel is formed and the color superconductivity, spontaneous breakdown of the color 
gauge symmetry SU{3)c SU{2)c, occurs This phenomenon can be understood 
in the sense that the effective dimensional reduction (3 + 1) (1 + 1) occurs with the 
large chemical potential as well as the case of QED4 in a constant magnetic field 0. 

There has been many studies of the color superconductivity in cold and dense QCD 
using the so-called Nambu-Jona-Lasinio (NJL) type models 0, the instanton liquid 
model 1^ and the Schwinger- Dyson approach to QCD in weak coupling limit 0, in other 
words, in yU cxo limit. The first two models have shown that the superconducting gap 
could be large as 100 MeV. On the other hand, due to the exchange of the magnetic 
gluons, the gap from QCD in weak coupling limit could be large as compared to that 
obtained under the condition which the magnetic gluons acquire static screening mass. 
If we extrapolate the gap down to the intermediate region of the chemical potential, 
it increases up to 100 MeV. However, as temperature grows, the color superconducting 



gap is reduced, then, the phase undergoes phase transition. The order of the transition 
seems to be second order in either case of with or without gauge fields. However, with 
gauge fields, there is a possibility of first order or a smooth crossover. 

We consider the color superconductivity under the situation where a flavor asymmetry 
takes place. The color superconductivity at finite flavor asymmetry has been studied in 
Ref . 0, . Especially, in Ref . , it has been argued that a mixed phase of the normal 
and the superconducting phase exists with a sufficiently large flavor asymmetry. In this 
paper, we study at finite temperature and discuss what extent does the mixed phase 
persist. Such an investigation has importance in studies of neutron stars and, of course, 
high-energy heavy-ion collisions at Alternating Gradient Synchrotron (AGS) and CERN 
Super Proton Synchrotron (SPS). 

The study of the color superconductivity in QCD is nontrivial and, therefore, we 
choose a tractable model to describe the quark-quark condensate. The most familiar 
of such models is known as the NJL (-type) model. Particularly, in the leading order 
of the 1/Nc expansion, the resultant gap function is independent of momentum and it 
makes our analysis simple one. Moreover, the structure of the color superconducting 
phase depends on the number of the quark fiavors. At extremely high density, where we 
can neglect the mass of s quark, the color-fiavor locked phase appears [§. However we 
expect, at relatively low density, there is a window of two-flavor color superconducting 
phase. Such being the case, we use the two-flavor, three-color NJL-type model. 

This paper is organized as follows. In Sec. II, we present our model and derive the 
thermodynamic potential at flnite density, temperature and flavor asymmetry. Then, in 
Sec. HI, we present numerical results for the thermodynamic potential. Section IV is 
devoted to the summary and discussion. Our unit are Ti = ks = c = 1. 



2 MODEL 

Our aim is to explore the phase diagram of two-flavor QCD. The SU{2)l x S'f/(2)/j 
chially symmetric Lagrangian [p!0[ 

^NJL = tpij ■ dtp + Cint, (1) 

Ant = Gi [(V5^)2 + (^^75^^)'] , (2) 

is our starting point. The quark fleld ip is assigned to the fundamental representation of 
the color SU{Nc) group and carries the flavor index i and the color index a as ipl^. The 
Pauli matrices r = {ti,T2,T3) act in the flavor space. In Eq. (|), we include two four- 
fermion interactions, namely, the scalar isosinglet (a) and pseudo-scalar isotriplet (vr). 
However, we do not abandon the possibility of adding other four-fermion interactions; 
e.g., the scalar isotriplet (oq) and the pseudo-scalar isosinglet (//'). 

Applying the Fierz transformations to these color singlet interactions, one can trans- 
form them into color 3 and 6 channels. In this paper, we consider the following scalar 
isoscalar color 3 diquark condensation which breaks SU{3)c to SU{2)c', i.e., taking a 
third direction preferred one in color space, 

A = As = a; ~ e^ha^sMfC^M, (3) 
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where C is a charge conjugation matrix, defined by C7^C^^ = — 7J and C = —C^^ = 
—C'^ and A is independent of momentum in our approximation (mean field approxima- 
tion or, equivalently, leading order in the 1/Nc expansion). The interaction Lagrangian 
which leads to the diquark condensation is written as 



dig 



Go 



(4) 



where is the antisymmetric Pauli matrix, is the antisymmetric (color 3) color 
generator, and the ellipsis denotes irrelevant diquark channels. Since the NJL model 
is not renormalizable, one cannot define the model until a regularization scheme has 
been specified. In this paper, we calculate the thermodynamic potential with a 3D 
sharp cutoff in momentum space. For a mesonic sector, we use a set of parameters: 
a coupling constant Gi = 5.01 GeV~^ and a 3D ultraviolet cutoff A = 0.65 GeV 
These parameters reproduce the value of = 93 MeV and the quark condensate (qq) = 
(— 225MeV)~^. The coupling constant G2, mediating an attractive four-fermion diquark 
interaction, is related to Gi through the Fierz transformations. In other words, the 
strength of G2 is fixed by the Fierz transformations. The Lagrangian (^ gives G2 = 
Gi/A. However, as mentioned before, one can change the ratio G2/G1 by adding other 
color singlet channel to Eq. (0). Therefore we take G2 independent on Gi and study the 
G2 dependence of A and a phase structure. 

We start our study with an investigation of the thermodynamic potential as a function 
of A at nonzero temperature, "averaged" chemical potential and "isospin" chemical 
potential. The averaged chemical potential /i^ and isospin chemical potential /xj are 
given as follows 



/is 



(5) 



where ^u{d) is a chemical potential of the u{d) quark. In order to evaluate the thermody- 
namic potential, one can use auxiliary field method with the saddle point approximation 
or the Cornwall- Jackiw-Tomboulis effective action |]ll| up to two-loop level with the 
Nambu-Gorkov formalism. These two methods are correspond to the evaluation at the 
leading order of the 1/Nc expansion and consistent with each other. The thermodynamic 
potential V{A] ^ib^I^i) in the leading order of the 1/Nc expansion is obtained as follows 



\/(A;/iB,/i/) 



A^ 

+ \n 
+2 In 
-F21n 
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(6) 



The finite temperature thermodynamic potential V{A;T, fiB, fJ^i) is obtained by using 
the formula: 

^^^.(.).ll/^/..tanh^.(.„), (T) 



(27r) 
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where (3 = 1/T and the a;-integration is to be understood that it is round an anti- 
clockwise contour including the pole of F{uj,q). Thus, we obtain the thermodynamic 
potential as follows 

+ Tin (1 + exp [-/3(e-(g) - fu)]) + Tin (1 + exp [-/5(e_(g) + /i/)]) 
+ Tin (1 + exp [-/3(e+(g) - fn)]) + T In (1 + exp [-/?(e+(g) + fn)]) } 

+ T In (1 + exp [~P{E4q) - fij)]) + T In (1 + exp + fn)]) 

+ Tin (1 + exp [-(3{EM " ^/)]) + Tin (1 + exp [-(3{E+{q) + ^,)])(^) 



where 



e±(g) = 1^1 ±/^B (9) 



are the gapless (color 3) quasi-particle energies relative to the (averaged) Fermi surface 
and 



E±{q) = {\q\±^lB), 



A2 

(10) 



\ (|q|±M^' 

are the gapful (color 1, 2) quasi-particle energies and 3D cutoff A is introduced to regu- 
larize the divergent integrals. 

The extremum condition for V{A] T, fiB, fJ^i) with respect to A lead to the Schwinger- 
Dyson equation (SDE) for A at nonzero temperature, density, and flavor asymmetry. The 
Cooper pairs condensate (ipip) is related Aq the solution of the SDE as 

iijij) = Ao/(2G2). (11) 

The "baryon number density" and the "isospin density" are obtained by 

d 

riB = nu + nd = — V{A;T, fis, fij) 

djiB 

d 

rii = nu-nd = — V{A;T, /xb, ^i)- (12) 

3 COLOR SUPERCONDUCTIVITY WITH ASYM- 
METRIC MATTER 

Figure 1 shows the G2 and fis dependence of A at zero temperature and zero flavor 
asymmetry. The value G2 = 3Gi/4 was taken by Berges and Rajagopal 0], where Gi is 
the coupling constant in a (scalar isosinglet) channel. Note that we use the Gi different 
from Berges and Rajagopal and 3D sharp cutoff instead of introducing a smooth form 
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factor F{q) = A^/(g^ + A^). Then, the behavior of A shghtly different from theirs. The 
value G2 = Gi/2 and G2 = Gi/A are introduced by the instanton vertex of two-flavor 
QCD and the SU{2)l x SU{2)r symmetric NJL model, respectively. As we can see, A 
is quite sensitive to (j2- For G2 = 3Gi/4, A could be as large as 100 MeV. However, it 
is about 50 MeV for G2 = Gxjl and negligible small for G2 = GijA. In the following, 
we use G2 = 3Gi/4 to get sizeable gap. 

Figure 2 shows the zero temperature and = 400 MeV thermodynamic potential as 
a function of A for several /i/. For A larger than /x/, all curves are coincident with each 
other. This is because, for A > /i/, the flavor asymmetry does not affect the thermody- 
namic potential. As /i/ grows, the potential in the trivial state (A = 0) decreases. It has 
two degenerate minima and the superconducting phase and the normal phase coexist at 
some value of /x/. For more large values of /x/, the normal phase is favored. As a result, A 
behaves constant up to a certain critical value of /i/ and, then, jumps to zero. We have a 
strong first order phase transition. We can see the same behavior of the thermodynamic 
potential and the order parameter in chiral symmetry restoration at zero temperature 
and finite baryon density. 

In Fig. 3, we present the phase diagram as a function of /is and /i/ at zero tempera- 
ture. Along the horizontal axis, color superconducting phase exists for arbitrary nonzero 

because the BCS instability exists for arbitrary small coupling strength. On the other 
hand, as mentioned before, we know that it undergoes a phase transition of first order 
with increasing flavor asymmetry. Then we get the phase boundary such as shown in 
Fig. 3. This boundary reflect the /i^ dependence of A at /x/ = (see Fig. 1). It is 
also interesting to see the phase diagram as a function of /is and /x/. Figure 4 shows the 
phase diagram in the Ub-uj plane. For the flavor asymmetry rij/uB larger than about 
15% the system remains in the mixed phase for arbitrarily high densities. This nature is 
almost the same as shown in Ref. 0. For the flavor asymmetry larger than about 40% 
the phase transition is complete and the system undergoes a phase transition into the 
normal phase. This point is differ from the results of Ref. and the difference arises 
from the choice of the parameters. 

We now turn to physics at finite temperature. Figure 5 shows the /// dependence of 
V for T = 20 MeV and Hb = 400 MeV. At relatively low temperature, a first order phase 
transition occurs. On the other hand, at relatively high temperature, the phase transition 
changes to into second order. From the behavior of the thermodynamic potential, the 
existence of a tricritical point in the T-fij plane is expected. In Fig. 7, we present the 
typical phase diagram as a function of T and fij for fiB fixed to 400/500 MeV. It is 
known that, in four-fermion interaction model which neglects the gauge fields, we have 
a second order phase transition for /x^ 7^ and /x/ = with increasing temperature 
(BCS theory predicts the critical temperature Tc is given by Tc ~ 0.57At=o)- On the 
other hand, we have a first order transition for T = with increasing flavor asymmetry 
/i/. Therefore the existence of a tricritical point in the T-fij plane , where the second 
order phase transition becomes that of the first order, is expected. In fact, we found 

T ~ 35 MeV, /i/ ~ 65 MeV for /i^ = 400 MeV, 
T ~ 42 MeV, fij ~ 69 MeV for fiB = 500 MeV, 
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from the condition 



a2(T,fij) = ai{T,fij) = 



(13) 



where we expand the thermodynamic potential, in the vicinity of the tricritical point, as 
a power series in A: 



Therefore, there is a possibihty of a first order phase transition as temperature grows, 
with sufficiently large flavor asymmetry. 



In summary, we studied the phase structure of the color superconductivity at flnite den- 
sity, temperature and flavor asymmetry by investigating the thermodynamic potential 
in the NJL model. 

Since specifying the ratio G2/G1 by the Fierz transformations in the mean-fleld ap- 
proximation is not unique, we examined the G2 dependence of A and found that A is 
quite sensitive to G2. 

First, let us see /i/ dependence of A. At zero temperature, as /i/ grows, the super- 
conducting gap A remains constant until fij is increased to its critical value and, then, 
vanishes. This behavior seems to be the same as obtained by Bedaque f^. Following 
Alford et ai, let us consider a concrete example; a gold nucleus at a baryon density eight 
times in nuclear matter. This matter, the number of d quark exceed that of u quark, 
corresponds to fiB = 535 MeV and = 12 MeV (up to sign). They have found that 
the gap is a little (about 18%) decreased. Our resultant gap does not changed in this 
condition. The distinction maybe arise from the difference in the coupling constant and 
the regularization scheme. 

Second, concerning the phase structure, we found the similar phase structure in the 
nB-ni/uB plane to that reported by Bedaque 0. Moreover, we found that the structure 
of the color superconducting phase is remarkably rich: 

(1) mixed phase appears at sufficiently large flavor asymmetry with temperature less 
than 40 MeV, 

(2) tricritical point exists in the T-/iB-/i/ phase diagram. 

The most favorable condition for the color superconductivity, cold and dense quark mat- 
ter, are at the core of neutron stars. However, the values of T, fiB and fij accomplished 
in heavy-ion collisions at, for example, AGS and SPS, maybe reach into the mixed phase, 
since it lies up until T ~ 40 MeV (see Fig. 7). Then, if the charged bubbles is formed, 
it may lead to observable effects, e.g., an obvious change of the vr"/??"*" ratio. 

Our results strongly depend on the regularization scheme, although qualitative fea- 
tures do not changes. In addition, it seems to be necessary to take into account next-to- 
leading order in the l/Nc expansion, because the ratio of G2/G1 has a order 0{1/Nc). 
Also, the dynamical chiral symmetry breaking and the pion condensation should be in- 
cluded in further study. The pion condensed phase has been conflrmed by the Monte 
Carlo simulation of two-color QCD |T^. We also plan to study the effective potential 



y (A; T, fxj) = \/(0, T, fxi) + a2{T, /x,) A^ + a,{T, m)A^ + ■■■. 



(14) 



4 SUMMARY AND DISCUSSION 
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within QCD in weak coupling limit and the three flavor case where the color flavor locked 
phase takes place. 
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Figure 1: Superconducting gap A as a function of /ib- The curves correspond to G2/G1 — 
3/4, 1/2, 1/4 (top to bottom). 
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Figure 2: The zero temperature thermodynamic potential as a function of A for several 
muj. The curves correspond to //j = 0, 50, 77.3, 100 MeV. The value fij = 77.3 MeV 
corresponds to a critical isospin chemical potential. Note that we neglect the contribution 
of color 3 quark to V which is independent of A. The same neglect will be done in Fig. 
5 and 6. 
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Figure 3: Phase diagram in the /ib-IJ'I plane. The (dotted) hne indicates the first order 
phase transition. 
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Figure 4: Phase diagram in the UB-rij/nB plane (in unit of lO^MeV^). 
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Figure 5: The thermodynamic potential for T = 20 MeV, /xs = 400 MeV and several 
Hi. The curves correspond to /// = 0, 71, 80 MeV (top to bottom). 
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Figure 6: The thermodynamic potential for T = 50 MeV, hb = 400 MeV and several 
Hi. The curves correspond to /x/ = 0, 49, 80 MeV (top to bottom). 
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Figure 7: Phase boundaries in the T-fij plane for fiB = 500 MeV (upper Une) and 400 
MeV (lower line). The solid lines describe the transition of second order and the dashed 
line describe that of first order. Points indicate tricritical points. 
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